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The question of the integrability of the mixmaster model of the Universe, presented as 
a dynamical system with finite degrees of freedom, is investigated in present paper. As 
far as the model belongs to the class of pseudo-Euclidean generalized Toda chains the 
method of getting Kovalevskaya exponents developed in P| for chains of Euclidean type, is 
used. The generalized formula of Adler and van Moerbeke for systems of an indefinite 
metric is obtained. There was shown that although by the formula we got integer values 
of Kovalevskaya exponents there were multivalued solutions, branched at particular points 
on a plane of complex time t. This class of solutions differs from ones considered in ||. 
Apparently, the system does not possess additional algebraic ||,|§] and one- valued first 
integrals because of complex and transcendental values of the exponents. In addition, in 
the section 4 a ten-dimensional mixmaster model |J was studied. There were not integer 
exponents. 



PACS numbers: 04.20, 04.40. 



1 The Pseudo-Euclidean Generalized Toda Chains 



The mixmaster model belonges to pseudo-Euclidean generalized Toda chains [p]] on a level 
H = 0. The hamiltonian has a form 

H = l -(-pl+pl+pl) + exp(4a)V(P + ,P_), (1.1) 

where the potential function V(/3+, /3_) is an exponential polinomial: 

V(J3+, P-) = exp{-8(3 + ) + exp{4(3 + + lyftpj) + exp(A(3 + - W3/3-) 

- 2exp(4f3 + ) - 2exp{-2(3 + + 2V3/3_) - 2exp(-2(3 + - 2V3/3_). (1.2) 
The hamiltonian has a form: 

1 N 
H = - < p,p > +Y, c i^> (1-3) 
1 i=i 

where <, > is a scalar product in a Minkowski space R l,n ~ 1 , C{ are some real coefficients, v% = 
exp(cii,q), (, ) is a scalar product in a Euclid space R n , a are real vectors. Pseudoeuclidity 
of a momentum space is a distinctive peculiarity of gravitational problems sothey can not 
be refered as analytical dynamics problems, where a quadratic by momenta form is a kinetic 
energy. 

On the other hand, the cosmological models can be considered as dynamical systems. 
So it is possible to carry over strict methods of analysis, traditionally used in the analytical 
mechanics, and adapt them to the systems like (1.3). For analysis of integrability of the 
cosmological model which is a hamiltonian system with three degrees of freedom, let us 
apply the Painleve test (see e.g. |7]]) for calculation of Kovalevskaya exponents ||. The term 
" Kovalevskaya exponents" was introduced in paper Q thus it was marked an outstanding 
contribution of the Russian woman into solution of the important problem of dynamics- 
integration of a task of a rigid body rotation. It was found that a general solution is 
meromorphic if an only if the system possesses an additional one-valued or algebraic [§] 
first integral. Further development of the theory of Kovalevskaya exponents is contained in 
0- 

By expanding 2n-dimensional phase space to 2iV-dimensional one by homeomorphism 
(p,q) h-> (v,u): 

Vi = exp(a h q), Ui=<a h p>, i = l,...,N, (1.4) 

generalizing the Flaska transformation [0 , one gets a hamiltonian system whose movement 
equations are system of autonomic homogeneous differential equations with polinomial right 
side: 

Vi = UiVi, 

N 

u l = Y,M ij v j , i = l,...,N. (1.5) 

3=1 

The matrix M is constructed of scalar products of vectors a in a Minkowski space R}' n ~ 1 : 

Mij = —Cj < ai,dj > . (1.6) 
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The property of integrability of dynamical systems appears in a character of singularities 
of solutions what it is not possible to say about points of common position, so just singular 
points represent particular interest for investigation. The equations (1.5) have the following 
partial solutions: 

Ui = Ui/t, v i = V l /t\ i = l,2,...,N, (1.7) 
coefficients Ui, Vi obey a system of algebraic equations 

-2Vi = UiVi, 

N 

- Ui = J2 M ijVj. (1.8) 

i=i 

Now let us analize all types of solutions of the system (1.8) in detail. 

I. Let V\ 7^ 0, the rest V 2 , V3, Vn = 0, then we get a solution: if Mn 7^ 0, then V\ = 
2/M 11 ,U 1 = -2,U 2 = -2M 21 /M n ,...,U N = -2M N1 /M n . Analogously the last solutions 
will be obtained. If for some % : Vi 7^ 0, and Vj = for all j 7^ i, then if Ma 7^ 0, and 
we get Ui = —2, Uj = —2Mji/Ma for all i ^ i. It follows from the obtained solutions that 
the significant point of analysis is a nonequality of a correponding diagonal element of the 
matrix M to zero that is possible in case of an isotropy of a vector a^. It is a principial 
distinctive feature of pseudoeuclidean chains. 

II. Let us pass to the next series of solutions. Let V\ ^ and V 2 7^ 0, the other Vi = 0, 
then U\ = U 2 = —2, V\ and V 2 are found from a system 

2 = MnVi + M 12 V 2 , 

2 = M 2X V X + M 22 V 2 . (1.9) 

The condition of existence of solutions of the system (1.9) is a nondegeneracy of a matrix 
M\ 2 . If detM\ 2 7^ 0, then, obtaining from (1.9) values V\ and V 2 and putting them in (1.8), 
one gets U 3 , U4, U N : 

-U 3 = M 31 V 1 + M 32 V 2 , 

- U 4 = M41K + M 42 V 2 , (1.10) 

-C/jv = M N1 Vi + MjvaVi. 

By the same way we obtain the rest of solutions of this series searching through all pos- 
sible nonzero pairs Vi 7^ 0, Vj 7^ 0. The corresponding condition of existence of solu- 
tions is a nondegeneracy of a two-dimensional matrix: detMij 7^ 0, constructed of elements 

Mjj, Mij, Mji, Mjj. 

III. The following series of solutions are obtained by considering all nonzero triple solutions 
Vi,Vj,Vk with a condition of nondegeneracy of corresponding three-dimensional matrices. 
The last possible solution is: U { = -2, i = 1, 2, N, V = —(M ij )~ 1 Uj, detM ^ 0. 
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2 A Calculation of Kovalevskaya Exponents 



For investigation of a single-valuedness of obtained solutions we use Lyapunov method [[| 
based on studying of behaviour of solutions of equations in variations. First variations obey 
the following differential equations: 

d N 

= — + -?-> » = L ■■■>*■ (2-1) 
We seek their solutions in the form of 

fcii = &^\ 5v i = r H t»- 2 , i = l,...,N. (2.2) 

Then for searching of coefficients 77^ one gets a linear homogeneous system of equations 
with a parameter p: 

(p - 2 - £7^ = 

iV 

(p-l)6 = £My»fc, i = l,...,N. (2.3) 

3=1 

Values of a parameter p are called Kovalevskaya exponents. 

I. Let us consider solutions of the first series when Vi 7^ 0. Let 771 7^ and the rest % = 0, 
then from the first system of equations (2.3) one gets £1 = Mnprji/2, a substitution of it to 
the second system (2.3) gives a condition on values of a parameter p:p(p — 1) — 2 = 0, i.e. 
pi = —1, P2 = 2. The remaining equation (2.3) gives us solutions £« = £1(771, p). 

Let then r/ 2 7^ 0, then 773,774, ...,^v = 0, p = 2 — 2M 2 i/Mn, pr/i = 2^/Afn. The second 
system gets functions £j = £4(7/2), i = 1, 2, /V. 

Further, considering 773 7^ 0, we receive p = 2 — 2 < a 3 , a! > / < a 1; a! >; r/ 2 , 774, r/ N = 
0; pr/i = 2£i/Mn; and from the second system: £j = £j(r/ 2 ), i = 1, 2, iV. 

As a result, having looked over all solutions of the first series for a case V\ 7^ 0, we obtain 
a formula for a spectrum p : p — 2 — 2 < a i: ai > / < a 1; >, i = 2, 3, iV. 

As a final result, having considered the rest solutions of the first series, we obtain a 
formula for Kovalevskaya exponents p that is the formula of Adler and van Moerbeke for the 
case of indefinite spaces: 

p = 2-2 <<li,ak> , i^k, <a k ,a k >^0. (2.4) 
< Ofc, a k > 

The requirement p E Z is & necessary condition of a meromorphity of solutions on a complex 
plane of t. It should be noticed that under derivation of the formula (2.4) no restriction on a 
metric signature was imposed. It is correct not only for spaces of the Minkowski signature. 

II. The following series of solutions do not leed to the similar elegant result. But we may 
show that solutions for variations exist in principle. So, for example, let us take a solution 
of the second series when V\ 7^ 0, V2 7^ 0. Then we have: 

PVi = Vi£i, 
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PV2 = V 2 &, 

(p - 2 - U 3 ) m = 0, (2.5) 



(p - 2 - U N )rj N = 0, 

N 

(p-l)& = E M ^> i = l,...,N. 

3=1 

Let us assume, for instance, 773 7^ 0, then 774, ...,t)n = 0, and therefore p = 2 + C/3. Out of 
the first pair of equations, expressing 77! = V"i£i/p, ?] 2 = V^/p, we substitute them to the 
second system. We obtain: 

(p(p - 1) - MnK)fi = M 12 \/ 2 6 + pM 13% , 

(p(p - 1) - M 21 K)^ = M 22 \/ 2 £ 2 + pM 23 r/ 3 , 
p(p - 1)& = M316 + M 32 £ 2 + Afepr/s, (2.6) 

p(p - 1)£ N = Mjvi£i + .Uv-jG + M N3 pr] 3 . 

The first two equations represent a system of nonhomogeneous linear equations, whence 
we get solutions £1 = £1(773), £ 2 = £2(773)- Substituting them to the rest equations, we get 
£3 = £3(7/3), £n = £^(773)- By the same way the other solutions are to be obtained. 

III. For the last solution when all Ui = —2 for % — 1,2, ...,N, and V* are solutions of a 
matrix equation 

N 

2 = £ MyV} (2.7) 
for the nondegenerate matrix M, we receive an eigenvalue problem: 

N 

EW-^-l))^ = 0- (2-8) 
j'=i 

As it is seen from (2.7), a sum of elements of any line of a matrix M^Vj is two. Out of 
this property of the matrix we obtain that one eigenvector £3 is unit with an eigenvalue 
p(p — 1) = 2. From here we obtain already known roots: p = —1, 2. 

3 The Mixmaster Model of the Universe 

Now let us apply the method to analysis of integrability of the mixmaster model of the 
Universe, "root vectors" of which have a form (see (1.2)): 

ai(4,-8,0), a 4 (4,4,0), 

a 2 (4,4,4V3), 05(4,-2,2^3), (3.1) 
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a 3 (4,4,-V3), a 6 (4,-2,-2>/3). 

"Cartan matrix", composed of scalar products of the vectors (3.1) in Minkowski space, then 
has a form: 



< a,- 



a,- >= 48 



/ 1 


-l 


-1 


-1 





\ 


-1 


l 


-1 








-1 


-1 


-l 


1 





-1 





-1 











-1/2 


-1/2 








-1 


-1/2 





-1/2 


V o 


-1 





-1/2 


-1/2 


/ 



(3.2) 



One gets, three "root vectors" dispose out of a light cone (space-like vectors), the rest 
three are isotropic, are situated on the light cone. Using the generalised Adler-van Moerbeke 
formula (2.4), taking account of zero norm of three vectors, we get p\ = 2, p 2 = 4, i.e. they 
are integer. 

Now we construct a matrix M by the formula (1.6), using a vector q(l, 1, 1, —2, —2, —2) : 



Mij = 48 



/-I 
1 
1 
1 


V o 



l 

-l 
l 




1 



1 
1 

-1 



1 





-2 




-1 
-1 






-2 
-1 


-1 



o\ 

-2 

-1 
-1 



(3.3) 



By virtue of a specific form of the matrix (3.3), it is not difficult a consideration of all 
partial solutions of the second type, when a pair Vi,Vj under some i,j is not zero. One 
may notice that there is only one, by virtue of nondegeneracy of blocks, nontrivial solution: 
V 1 ^ 0, V 4 ^ 0. It has a form: VI = 1/24, V 4 = -1/24, \J i = -2, i = 1, 6. 

Equations on coefficients of variations of solutions are: 

(p(p-l) + 2)77! + 4774 = 0, 
27/! +p(p~ 1)^4 = 0, 

(3.4) 



(p- 


1)6 


= M 2 i77i, 


(p- 


1)6 


= M 31 rj 1: 


(p- 


1)6 


= M bA 7] A , 


(p- 


1)6 





The condition of existence of solutions of the system of the first pair of equations is a 



factorized algebraic equation of a fourth order on to a spectral parameter p : 

[p(p-l)-2][p(p-l)+4]=0. 



(3.5) 



The solutions of the first square equation are integer: p 1 — — 1, p 2 — 2, but roots of the 
second equation 

p li2 = (l±i>/l5)/2 (3.6) 
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are complex and irrational. Thus, in a common case, a solution is not one-valued on the 
complex plane of time t. Under some additional conditions |§ , [f| , || , it leads to nonexistence 
additional algebraic or one-valued first integrals. Let us remark that meromorphity of a 
common solution is a criterion of existence of an algebraic first integral in dynamics of a 
rigid body (Gussont theorem) and in a problem of three bodies (Brunce theorem). 



4 A Mixmaster Model with Geometry R 1 x S 3 x S 3 x S 3 

In paper [fj a model of a spatially homogeneous and isotropic vacuum universe with geometry 
R 1 x S 3 x S 3 x S 3 was studied. A ten-dimensional model was chosen due to promising 
advances in a ten-dimensional superstring theory. Its hamiltonian in Misner coordinates 
(a,P±,9,(j>±,T},il;±) is: 

H = -f| + £ f + e -^[e-' e/V ~ 3 9W) + e 2 ^(e~ 2 M0) + e'MM, (4-1) 
where g(x) means a function 

g( x \ = e 4x++4,V3x- e 4x + -iV3x- e ~8x+ _ 2 e 4x + — 2 e 2x ++ 2 ^ x - — 2 e 2x +~ 2 ^ x - . (4.2) 

The problem also belongs to the generalized pseudo-Euclidean Toda chains. Its potential 
energy is a sum of 18 terms: 

18 

U = ^2ciexp(a h q). (4.3) 
i=i 

So 18 "root vectors" a in 9-dimensional Euclidean space are: 

ai(16, 4, 4^3, -4/V3, 0, 0, 0, 0, 0), 

a 2 (16, 4, -4V3, -4/V3, 0, 0, 0, 0, 0), 
a 3 (16, -8, 0,-4/^3, 0,0, 0,0,0), 
a 4 (16,4,0,-4/V3,0,0,0,0,0), 
a 5 (16, 2, 2^3, -4/V3, 0, 0, 0, 0, 0), 
a 6 (16, 2, -2^3, -4/V3, 0, 0, 0, 0, 0), 
a 7 (16, 0, 0, 4, 4^3, -2, 0, 0), 

a 8 (16, 0, 0, 2/V3, 4, -4V3, -2, 0, 0), 
a 9 (16, 0, 0, 2/y/S, -8, 0, -2, 0, 0), 
a 10 (16, 0,0, 2/^3,4,0,-2, 0,0), 
an(16, 0, 0, 2/y/3, 2, 2V3, -2, 0, 0), 
012(16, 0, 0, 2/y/3, 2, -2y/3, -2, 0, 0), 
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ai 3 (16, 0,0, 2/^3,0,0,2,4,4^), 
a 14 (16, 0, 0, 2/y/S, 0, 0, 2, 4, -4^3), 

ai 5 (16,0,0,2/V3,0,0,2,-8,0), 
a 16 (16, 0,0, 2/^3,0, 0,2, 4,0), 

ai 7 (16, 0,0, 2/^3,0,0,2,2,2^3), 

a 18 (16, 0, 0, 2/V3, 0, 0, 2, 2, -2^3). 

As before, we pass to redundant coordinates. Then get (l+8)-pseudo-Euclidean space 
with metric 77^ = diag(— 1/12, 1, 1). By the same way as in section 3 we calculate Ko- 
valevskaya exponents using the Adler-van Moerbeke formula (2.4). But besides integer 
numbers p = 0,2,3,4 there are rational exponents p\ = 4/3, P2 = 10/3. Apparently the 
considered model is nonintegrable. 



5 Discussion of Results 

A search of first integrals in the form of algebraic and one-valued functions is too restrictive. 
A more interesting question is about existence of additional real-analytic integrals. Such a 
statement of the question goes back to Poincare. For Euclidean Toda chains a question about 
existence of additional real-analytic polinomials by momenta first integrals was studied in 



paper |H| where a classification of integrable systems was brought. 

Furthermore, in this case, we are not interested in integrability at the all phase space but 
at a concrete hypersurface "integral of energy" H = 0. 

The next step of investigation of the presence of chaotic characteristics is a consideration 
of "truncated" model which to be obtain if one of the degrees of freedom of the mixmaster 
model is fixed. I.e. it is obtained if we rewrite a potential part of the hamiltonian (1.2), 
introducing coefficients in front of its terms 

V(P + , P-) = C ie xp(-8(3 + ) + C 2 exp{A(3 + + W3/3_) + C 3 exp(4p + - A^(3-)- 

2C i exp{A(3 + ) - 2C 5 exp(-2(3 + + 2>/3/3-) - 2C 6 exp(-2(3 + - 2V3/L), (5.1) 

and put C2 = C3 = C5 = C*6 = 0. Such a hamiltonian system in a four-dimensional phase 
system can be studied with the help of the numerical analysis of Poincare mapping if its 
phase trajectories are situated in a compact region. An existence of chaos in the "truncated" 
model implies a stochasticity of the initial model. Let us remark that a calculation by the 
generalized formula of Adler and van Moerbeke (2.4) for this case yields integer p, as to 
additional exponents, they are also complex and irrational, coincide with (3.6). 

Because the Killing metric in the generalized Adler-van Moerbeke formula is indefinite it 
should be pointed out at a classification scheme of noncompact Lie algebras by analogy with 
using the Cartan scheme of classification of compact simple Lie algebras for getting exact 
solutions of Toda lattices as it was done in |L4| . 



It is worth to mention paper [13], where was made an attempt of proof of complicated 
behaviour of the mixmaster model trajectories out of an analysis of level lines of the potential 
energy function. 
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